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ABSTRACT 

Inverse cascades of magnetic quantities for turbulent incompressible 
magnetohydrodynamics are reviewed, for two and three dimensions. The theory is 
extended to the Strauss equations, a description intermediate between two and 
three dimensions appropriate to tokamak magnetofluids. Consideration of the 
absolute equilibrium Gibbs ensemble for the system leads to a prediction of an 
inverse cascade of magnetic helicity, which may manifest itself as a major 
disruption. An agenda for computational investigation of this conjecture is 
propos'id. 
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A, Background 


I . INTOODUCTION 


Over the last fifteen years, the term "inverse cascade" has been ap- 
plied to a class of turbulent processes in which the small spatial scales 
Ihigh waveniunbers ) feed some q.uantity, by nonlinear processes, to the large 
jpatial scales (low wavenumbers). A variety of quantities can be inversely 
cascaded, depending upon the medium and upon the geometry; the process can look 
iuite different in different cases. Inverse cascades have been most intensively 
studied in Navier-Stokes fluids and incompressible magneto fluids. Most progress 
has been numerical, and the accurate computation of an inverse cascade is a 
strenuous numerical chsillenge because of the high spatial resolution required. 
This needed high spatial resolution considerably restricts the range of allowable 
Reynolds numbers , and the implications of these restrictions have yet to be 
felt in the plasma simulation commimity, though the hydrodynamicists are pain- 
fully aware of them. 

The idea of an inverse cascade first appeared in connection with flows 


in two-dimensional' 




Navier-Stokes fluids, then in nagnetofluids , first 


in three dimensions^ *"^*®*^ and then in two.^*^’^^*^^ Most of the work has assumed 
the simplest and most highly idealized boundary conditions: rectangular periodic 

ones. But now, possible relevance to the "magnetic dynamo" and to magnetic 
fusion confinement is exerting pressure for computations with fewer idealiza- 
tions and more realism. Computations sure crucial for demonstrating most of the 
theoretical consequences. Analytical indicators can be used to anticipate the 
existence of inverse cascades, but accurate prediction of any of their details 
seems inherently to require sophisticated numerics. No one has had much success 
proceeding solely with pencil and paper. The intention here is to clarify the 
implications of certain classes of possible computations. 
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The mechanism underlying an inverse cascade is essentially simple. * 

Two things seem to be necessary: (l) a dissipation which only becomes effective 

at small spatial scales, such as viscous or Ohmic dissipation; and (2) two or 
more non-dissipative integral invariants of the motion which are representable 
us sums which emphasize differently the large and small scales. 
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Take an example: tvo-dimensional Navier-Stokes flow at higti Reynolds 

numbers (the same mathematics does as well for electrostatic guiding center 
plasmas). The two non-dissipative invauriants, expressed in terms of the Fourier 
coefficients of the velocity field y(k,t), are the energy per unit mass E and 
the mean-square vorticity or "enstrophy” f2: 

The non-dissipative terms in the equations of motion are the nonlinear ones, 
and if they transfer excitations to the high ]k| part of the spectinan (as they 

surely must), the only way the simultaneous constancy of both E and can be 

preserved is for there to be an accompanying transfer to low jk|. The additional 
constant of the non-dissipative motion Q (which is not conserved in three di- 
mensions) demands that the non-dissipative (i.e., non-linear) terms in the 
Navier-Stokes equation must transfer excitations in both directions in wavenumber 
space, if there is any transfer at all. The paths in k space by which the transfer 

occurs are a secondary issue, but there is nothing simple about them. 

Subtleties arise when one begins to ask sharper questions, such as 
which quantity gets transferred in which direction and how fast. These questions 
have been earnestly addressed in the literature and do not lead to short, graphic 
answers, only to long, contingent ones. Relatively simple and persuasive answers 
can be given for the quasi-steady state, when a source of excitations, band- 
limited in wavenumber space, is regarded as injecting the excitations at a sta- 
tistically steady rate and attention is directed toward the spectra and the 
transfer rates of the cascaded quantities. The word "source" can mean lots of 
things, from an externally-applied forcing field (such as an electric field or 
current field, in magneto-hydrodynamics) to some microscopic instability which 
might be present. 

Computationally, one often works initial value problems starting from 
smooth initial conditions, and in that case, it may be useful to identify a 
particular unstable normal mode as associated with the onset of the turbulent 
fluctuations. In the laboratory one usually does not work an initial value 
problem starting from smooth initial conditions; rather, unstable systems are 
often created only with some level of excitations already in place. It may then 
be a less interesting question as to which "mode" is the "most unstable" one. 
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because the fully developed and observed state may bear little relation to the 
kind of smooth profiles on which any kind of stability calculations can be done. 

A hydrodynamic analogy is pipe flow far above the critical Reynolds number; 
there is hardly even an academic relation between the details of the flow and 
any stability calculation that can be done. Particularly in plasma confinement 
experiments, there may be many simultaneous sources of excitation: a gradient 

in virtually any mean field variable (pressure, magnetic field, electric current, 
flow velocity, density, temperature) is potentially a source of excitations. 

Many sources may act simultaneo\isly. They may interfere with or reinforce each 
other. 

What happens to the excitations once they are launched on their journey 
through k-space may be the more interesting and the more physically significant 
question. There is a clear difference between turbulent processes that transfer 
their distxurbances to small rather than large scales from the point of view of 
thermonuclear confinement. In the former case, the worst that can happen is an 
enhanced transport, a disadvantage which may well be offset by an enhanced heating. 
(A tokamak which had no ’’directly cascaded" turbulence of this kind might not heat 
at all! ) Transfer toward large scales can clearly be more serious: macroscopic 

bulk deformations can redi.stribute and disrupt the plasma variables and perhaps 
terminate the confinement. Disruptions in tokamaks have some of the flavor of 
what one might expect from an inverse cascade, though experiments with highly 
resolved enooigh diagnostics to study the small scales of what happens in a major 
iisruption appear not to have been carried out. (A respectable beginning on the 
pr“blem was made some time ago by Hutchinson^^ and Morton.^^) For the near-term 
future, connections between tokamak disruptions and possible inverse cascades are 
likely to continue to be primarily a nximerical subject. 



B. Numerical Examples 


We briefly remark upon three examples of attempted computations of 
inverse magnetic cascades. The quantities expected to be inversely cascaded 
are mean square vector potential in two dimensions and magnetic helicity in 
three. 
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Figrxre 1, due to Pouquet, displays the results of a closure calcula- 
tion of an inversely cascading magnetic vector potential spectrum in two dimen- 
sions. The mechanical and magnetic Reynolds numbers are of the order of ii500. 

The excitations are injected at a wave number indicated by an arrow, and every- 
thing at lower wave numbers is the result of back transfer in wave number space. 

The computation is not a direct solution of the dynamical equations themselves, 
but rather a statistical, eddy-damped "closure” approximation to them which per- 
mits higher values of the Reynolds numbers than a direct computation could possibly 
permit with present-day computers. 
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Figure 2, due to Fyfe et al., shows a direct numerical solution of 

the two-dimensional magnetohydrodynamic equations, driven by a random magnetic 

2 

forcing confined to the indicated wavenumber band. All the points below k = 55 
correspond to inverse magnetic transfer to long wavelengths. This, incidentally, 
is a calculation in which spatial resolution was inadequate to the Reynolds number 
chosen, and as a consequence there is no "dissipation range", or region of pre- 
cipitous fall-off at increasing k, apparent at the upper end. 

IT 

Figure 3 is tedten from Meneguzzi et ail., and illustrates the results 
of a direct con^iutation from the three-dimensional incompressible magnetohydro- 
dynamic equations, again with a random forcing, but this time with a helical., 
mechanical one. A weak "seed" magnetic field amplifies with time as the injected 
kinetic energy converts into magnetic energy through dynamo action. The back 
transfer of magnetic helicity is apparent. 

Additional direct inverse cascade computations which may be mentioned 

are those of Lilly^^ and Fyfe et al.,^^ for the two-dimensional Navier Stokes 
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case; see also Pouquet et ail. for the corresponding closure computation. Pouquet 

7 

and Patterson have displayed three-dimens ionaLL magnetohydrodynamic direct compu- 

g 

tations. Pouquet et al. have given three-dimensional inverse cascade closure 
computations for incompressible magnetohydrodynamics. 

All the computations cited assume rectangular periodic boundary conditions 
with no net flux of any quauatity through any cross section of the system. 



II. THE STRAUSS EQUATIONS; ABSOLUTE EQUILIBRIUM ENSEMBLE THEORY 


Rectangular periodic boundary conditions, which characterize most of 
the computational efforts to study inverse cascades, idealize away several key 
features of real systems. In particular, they preclude net fluxes of such quan- 
tities as electric current through a cross-secticn of the region of computation. 

They also rule out all true boundary effects, which undoubtedly play a role in 
real situations. It is desirable to state the ideas associated with inverse 
cascades in a mathematics that corresponds to more realistic representations of 
actual plasmas. 

20 

The Strauss equations are a set of magnetohydrodynamic equations that 
are close to those of incompressible two-dimensional magnetohydrodynamics, but 
which include some important three-dimensional effects. They are far more tractable 
than the full set of three-dimensional magnetohydro dynamic equations. What ap- 
pears to the writer to be a more transparent derivation than Strauss's can be 
given, ending with the same equations; the derivation is not simple, however. 

Their most important feature is that of their near two-dimensionality: the 

variable magnetic fields and velocity fields are perpendicular to the z-direction 
(say), but the field variables eu:e functions of all three spatial coordinates. 

The variable part of the magnetic field is B and is expressed in teims of a vec- 
tor potential A as B = x (The subscript "A” will always mean perpendi- 

cular to e^.) The velocity field v is expressed in terms of a stream function 
U as V = Vj_ X © 2 ^* There is a constant, uniform dc magnetic field in the z 
dii 'ction, of magnitude Bp. » |b 1 . The vorticity to and the ctirrent density j are 

in the z direction and are given by V, U = -to and V, A = - j , in terms of A and U. 
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The Strauss equatio.ns are, in a familiar set of dimensionless variables 
in which flow velocities sure measured in units of the Alfven speed, 

II ♦ v.7,A - = B, f (1) 

p (It * 1-^“) - s- V - H <=> 

The uniform dimensionless mass density p can be consistently set equal to unity. 

The quantities and v ^ are essentially the magnetic and mecheinical Reynolds 
numbers, respectively. 
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Except for the terms on the right hand sides of Eqs. (1) and (2) and 

the z-dependences , Eqs. (l) and (2) are identical with the equations of two- 

dimensional incompressible magnetohydrodynamics, a system which is now 

gecting to be rather well understood. Many Interesting questions can be asked 

of Eqs. (1) and (2). For example, the question of possible stability thresholds 

for quiescent, slowly-decaying equilibria as y and \) decrease towards zero is 

an analogue of the question of the stability of hydrodynamic pipe flow or Couette 

flow as the Reynolds number increases. Computational studies of Eqs. (l) and (2) 

21 22 23 2k 

have been initiated by Carreras, Hicks and others » » » in connection with 

tokamak confinement. 

The following considerations are intended to facili' ite this enterprise 

by sharpening the mathematical framework in w..ich the interpretation of a major 

disruption as an inverse magnetic cascade may be niunerically tested. In a paper 

12 

published in 1977* put forward the suggestion that the inverse cascade of 

magnetic helicity might be responsible for major tokamaJc disruptions, and a test 

of this hypothesis seems now to be within re^'ch. Related considerations are 

25 

discussed in a forthcoming paper by Tetreault. 

All inverse cascade computations to date have been motivated by con- 
sidering first a model problem in which the non-dissipative dynamicaJ. equations 
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have been expanded in a set of orthogonal functions. ’ The expansion is then 
truncated at a large but finite number of terms. Statistical mechanical proce- 
dxires are then performed in the phase space defined by the expansion coefficients. 
Somewhat surprisingly, absolute equilibrium canonical ensembles (Gibbs distribu- 
tions) have proved to be accurate predictors of time averages of phase functions 
of the expansion coefficients. Considering the limiting behavior, as the number 
of expansion coefficients becomes infinite* highlights any tendencies which may 
exist for some invariant to migrate to long wavelengths. 

The Gibbs distributions are constructed from those non-dissipative 
invariants which remain invariant after the truncation: i.e., are "rugged”. For 

the Strauss equations there are apparently three such inigged invariants. They 

are the energy E, the "cross helicity" H , and the magnetic helicity H : 

c m 

E = /rdrdBdz + y^] (3) 

H = /rdrdSdz caA (U) 

c 

H = B^/rdrd9dz A. (5) 

m u 

The integrations are over the three dimensionsil region 0^9 < 2ir, 0j<r < a. 
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0 < z < L , where periodicity in z with period L is assumed. Free~slip boundary 
^ z z 

conditions vS = 0, B*§ =0 are assumed at r = a, the wall of the rigid cylinder, 

r ^ r 

For added realism (but added complexity), we might, add the "no slip" boundary 
conditions on the tangential components : x B) x = 0 and v x = 0. These 

are superfluous, however, for the model problem at hand, and greatly complicate 
the mathematics. 

Usually, only quadratic invariants stre easily proved to be "rugged", 
and this can be achieved by the somewhat formal device of treating B^ itself as a 
phase space coordinate whose equation of motion just happens to be dB^/dt = 0. 

Bq itself is then a rugged invariant, and all the terms in Eqs. (l) and (2) as 
well as the expressions (3)-(5) can be considered as quadratic. The ensemble 
chosen can be chosen to be sharp in Bq! i.e., microcanonical in Bq but canonical 
in the other invariants . 

A and U are expanded in the complete orthonormal set of eigenfunctions 
of the Laplacian,^^ 

nmq nmq m nmq n 

where k = 27rn/L , m and n ajre integers, and Y a is the qth zero of J (x). 
n z * 0*1 jjjjjq -1 jjj 

The normalization constant C is 

nmq 

C = (ttL a^)"^/^/J (Y a). (7) 

nmq z m+l nmq 


If we define + k V^A = -X^ A , and V.^A = -y^ A 

rimq nmq n 2 ®®*! 2 nm^l “mQ. 

■^nmq eigenfunctions of both and V . 

We may write the infinite sums 


u 

= Z 

n A 


nmq 

nmq nmq 

A 

= Z 

5 A 


nmq 

nmq nmq 

E 

= Z 

Y^ ( 1C_ 


nmq 

'nmq ' nmq ' 

[ 

= Z 

y" 5* 

c 

nmq 

nmq nmq nmq 


H = 2B 
m < 


0 (tf 


' nmq ' 

^ / complex V 
^ conjugate^ 


^=AZ^ 


where for economy of notation, A = 2BQ(7r/L^)‘^' The sxmis are over large but 

finite sets of terms once the truncation is performed. Also to be truncated is 

the set of ordinary first-order differential equations for dC /dt and dri /dt 

nmq nmq 

that result when Eqs. (8) and (9) are inserted in Eqs. (l) and (2), 


'OOq 

1/2 
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The Gibbs ensemble which is appropriate to the case of no initial cor- 
relation between v and B is the multivariate probability distribution 

= const. X exp {-oE - 8H } , (13) 

eq. ^ m 

with reciprocal temperatures chosen to match desired ensemble expecta- 

tions <E>, <H >. 

’ m 


Inserting Eqs. (10) and (12) into (13), the modal expectations are 
readily calculated. First 


p 2 

>=<n >=0,m +n ^0, and 

nmq nmq * ’ 


<| 5 , 


Also 


nmq' 


■> = < h 


■> = 


*^^00q^ 2a 


nmq' 

SA 1 


m^ + n^ 5^ 0 . 




nmq 


'OOq 


- <5ooq’>"’ - 


(lu) 

(15) 

(16) 


“^OOq 

a and 3 are determined as the roots of 


<E> = 


2 i + 2 + 


nmq 

n^+m^j^O 


a 


q a 


teV 

’^OOq 


1 


and 


m 2a 


l~r- 

^OOq 


(17) 


(18) 


which keep all the <|C I > and <|ri I > positive. 

nmq nmq 

If <E> and <H > are held fixed and the number of nmq modes is allowed 
m 

to increase without limi'., it is easy to show that a l3l with |a/3l a 

finite ratio. The sum in Eq. (l8) is convergent as the maximum q Except 

2 2 

for the excess energy distributed in infinitessimal increments over the m + n #0 


modes, the fl.id excitations freeze into the OOq magnetic modes, 
potential <A> approaches a function 


The vector 


<A> -*■ Z <E > A 

q=l ^OOq OOq 

_ BA 2. -1/2 ” 

= - 2^ ^ qil 


^OOq'^l^^OOq®^ 


a universal function whose defining expression is manifestly convergent, 

i2. 


(19) 


For 


all m^ + n^ 0, 


<K, 


nmq ' 


<£ >‘ 
^OOq 


as the number of terms approaches infinity. 


This behavior is slightly different from the corresponding inverse cas- 
cade behaviors found in previously examined cases, insofar as the helicity does 
not condense into the 001 mode alone but into the rapidly converging series de- 
fined by Eq. (19). The state (19) is in fact a uniform cxirrent density state , 

2 

with <A> varying proportionately to r . Nevertheless, the earmarks are there, 
and to anyone familiar with the previous history of the theory of inverse cas- 
cades, provide a basis for conjecturing that in the presence of dissipation and 

forcing, there would be an inverse cascade of H to long wavelengths, headed 

m 

for the state defined by Eq. (19). 

The uniform-current state (19) has the additional significance of being 

either the state of minimum energy for given helicity, or the state of maximal 

helicity for given energy. This can readily be seen from the Euler equation for 

2 

the variational problem of minimizing E subject to a fixed value of H ; V,A = 

m * 

a constant. For the dissipative initial value problem, it has been previously 

demonstrated that under many circumstances, the directly cascadable invariants 

will "selectively decay" relative to the inversely cascadable ones, and their 
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ratios will approach their theoretical lower bounds for large times. * For 

the Strauss equations, it appears that this loniform current state is the 
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"selectively decayed" state, analogous to the "Taylor state" or force-free 
state of the full set of magnetohydrodynamic equations. 

If the more realistic "no slip" boundary conditions are invoked, re- 
quiring J = 0atr = a, the Jiniform current state is not attainable. At the 
least, a sharp current gradient must develop somewhere across the cross-section 
of the cylinder as the selective decay progresses, perhaps in the form of a 
boundary layer near .r = a. (Something similar was seen in a recent selective 
decay calculation inside a compact toroid. ) This current gradient, which 

necessarily exists at the edge of a current-carrying plasma boimded by a con- 
ducting wall, looms as a rather universal and difficult-to-avoid source of 
"tearing mode turbulence." 

This gradient in the current near the walls is one among many poter- 

tial soiurces for the helicity which might be inversely cascaded. There is 

probably no single mechanism for supplying small-scale helicity. A second 

12 

likely possibility is current filament at ion that may develop along loced. hot 


spots in the m agnetofluid ; the resistivity falls off with increasing temperature. 
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and may thus channel the current alorg hot tubes of force, resulting in still 
higher local heating. There is probably no single mechanism for supplying 
smai,l-scale helicicy, which basically results any time current flows along a 
field line in the presence of resistivity. Arguing in favor of particular 
drivers or sources may prove as fruitless an activity as the generation of 
debates that has surrounded linear instability theory. It appears imperative 
to produce major disruptions in the presence of as many parameter variations 
as possible, to begin to acquire discrimination among the variety of sources 
that may be operative. 



III. A POSSIBLE COMPUTATIONAL AGENDA 


Numerical demonstration of the possible inverse cascade properties of 
Strauss's equations will be facilitated by taking advantage of earlier experience 
and conceptualizations gained in studying inverse cascades, and by resistance to 
letting the conceptual framework be circumscribed by linear stability analysis. 

The first limitation that will have to be confronted will be the limitation on 
the Reynolds numbers A currently popu3.ar rule of thumb for fluid com- 

putations is that for every unit of Reynolds number, one grid point (or finite 
element, or expansion coefficient) is required in each spatial dimension. Thus 

O 

a three-dimensional simulation at a Reynolds number of 30 requires about (30) 
grid points to resolve the smallest spatial scales. This limitation might be 
violated by a factor of two, but probably not by an order of magnitude. At 

present the smallest p and v that are feasible to compute with are between about 

-2 -3 -1 -2 

10 and 10 for two dimensions and 10 to 10 for three. The desired physical 

values, for real experiments, are likely to be considerably smaller than that 

(particularly p) . !Chere is no simple way around this difficulty. 

Progress can be made in perhaps only one of two ways. (1) It may be 

attempted to escape the connection, between the aforementioned necessary 

spatial resolution and wave number requirements by introducirg an artificial 
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enhanced dissipation which only becomes effective at high wave numbers. Ihie 
hope is that small scale dissipation only provides a sink anyway, and that the 
large scale dynamics will become independent of the details of the sink. (2) Al- 
ternatively, one may settle for the qualitative demonstration of the physics and 
extrapolate crudely to the behavior at the very high R'^ynolds numbers. In the 
former choice, the ultimate validation can only be a comparison with the results 
of very high-resolution codes in which no anomalous dissipation is introduced. 

In the second choice, the ultimate recoiirse is most probably to experiments. 

The most single important possibility, as far as major disruptions are 
concerned, has to do with the possible existence of thresholds in Reynolds num- 
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bers. Such thresholds, at sometimes surprisingly low values, have characterized 
the three-dimensional, computations on the dynamo problem: above certain critical 

Reynolds numbers, a given stirring mechanism will initiate an inverse cascade. 
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otherwise it will not. Of particular interest to tokamaks is the possibility 
that because of the drop of resistivity with increasing temperattire , critical 
Limdquist numbers (Alfven speed x length sc ale /magnetic diffusivity) or magnetic 
Reynolds numbers ( flow speed x length scale/magnetic diffusivity) will be crossed 
as the magnetofluid heats up, and an inverse cascade will start dramatically. 

A spectral -method numerical solution in search of inverse heli'sity 
cascades might proceed as follows, for Eqs. (l) and (2). Pick u, v 'v 1/50 or 
so, with a resolution of the order of, say, 100 in the radial and azimuthal 
coordinates and perhaps 15 to 20 in the z coordinate. (It is unknown as to how 
rapidly small scales in z will multiply themselves, but if they develop as 
readily as the small scales in 9 and r, the Strauss equations probably are not 
useful, anyway.) There is no reason to believe less spatial resolution is re- 
quired in 0 than in r. Invariance of the resiilts to increased resolution in 0, 
r, and z is a necessary check on their accuracy. An initial current profile whi -h 
is thought to be close to experimental reality should be chosen and an external 
forcing term (probably from a random number generator) should be permitted to 
drive the magnetofluid at the small scales, either locally or randomly in space 

as well. This can be most effectively accomplished by adding a small random 

l8 

term to the right hand sides of Eqs. (l) and (2), in the manner of Lilly, or 
12 

Fyfe ^ al . 

The purpose of this first exercise would be simply to see if an inverse 
cascade with the featiires of a major disruption can be induced, either by lowering 
the dissipation coefficients or raising the strength of the random forcing. Once 
such an event has been shown to exist, an infinite variety of refinements of the 
calculation, such as eliminating the random forcing in favor of a resistivity 
which depends upon the local temperature T, are imaginable. T, for example, 
might be convected and grow locally due to Ohmic dissipation: 

T = kv/t + nj^, 

where n is a heating rate, and K is a thermal conductivity. 

Finally, the most important question connected with major disruptions 
is not so much can they occur as what can be done about them. If their inter- 
pretation as an inverse cascade of helicity is correct, then at a formal level. 
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the remedy is apparent; feed the magnetoflui,d some helicity of opposite sign. 
This means in effect inducing in the plasma some current flowing in the opposite 
direction from the main toroidal current. Where and how to do this seems like 
a delicate matter. Not much help is to be expected from attempts at suppression 
by static external helical windings, which supply a vacuum helicity at a parti- 
cular n,m mode, but do little to the bulk distribution over modes that provides 
the basis for an inverse helicity cascade. 

The clear-cut demonstration of an example of inverse cascade behavior 
for the Strauss equations would londoubtedly stimulate many additional refinements 
and insights that are now hard to foresee. The time is overdue, also, when the 
plasma simulation community should begin to seek a conceptual framework for its 
computations in their hydrodynamic antecedents, and cease to expect that linear 
stability analyses will provide adequate insight into processes in which non- 
linear transfer is the dominant effect. 
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FIGUEE CAPTIONS 

Fig. 1 . Magnetic potential spectrum as a function of wave number for successive 
times, in two dimensions. Magnetic fluctuations and kinetic energy are injected 
at the wave number indicated by the arrow, and their transfer is approximated by 

•j ^ 

an eddy-damped "closure" calculation. (Taken from. Pouquet. ) 

Fig. 2 . Magnetic modal energies, averaged over angle, as function of |k|, for 

a direct solution of the two-dimensional magnetohydrodynamic equations in two 

12 

dimensions (taken from Fyfe et al. ). Magnetic fluctuations are randomly 
injected in the band between the arrows, and the spectrum is allowed to fill 
up. Initially, it is empty. 

Fig. 3 . Omni-direction magnetic energy spectrum for the three-dimensional case 

17 

(taken from Meneguzzi et al. ). A small seed magnetic field is amplified by 

2 

a driven velocity field which contains mechanical helicity, /y*0) d x. The 
spectra are labeled by the appropriate values of the time. This is the most 
clear-cut computation of "dynamo" action to date. 
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